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We observe the formation of momentum distributions indicative of spatial period-doubling of su- 
perfluid Bose-Einstein condensates in periodically translated optical lattices. The effect is attributed 
to dynamic instability of the condensate wavefunction caused by modulation-induced coupling of 
ground and excited bands. 
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Atomic Bose-Einstein condensates in optical lattice po- 
tentials have been a focus of many studies in the last 
decade, serving as a model system for study of tradition- 
ally solid state phenomena. Interparticle interaction, in 
combination with the periodicity of the lattice potential, 
give rise to interesting nonlinear dynamics, such as gap 
soliton formation Q and instabilities Hi These have 
been successfully described by the Gross-Pitaevskii equa- 
tion (GPE), which describes the mean-field evolution of 
the condensate wavefunction. Stationary solutions of the 
GPE have been shown to include Bloch waves, as well as 
non-Bloch states with longer spatial periodicity 0, and 
have been predicted to be energetically and dynamically 
unstable for specific quasimomenta ii Instability oc- 
curs when the excitation spectrum of the condensate is 
altered by the periodic potential and interaction in such 
a way as to make the macroscopic wavefunction unstable 
to perturbation. This has been experimentally demon- 
strated with condensates accelerated relative to the lat- 
tice @,0- 

We report here on the observation of a period-doubling 
dynamic instability introduced by periodically translat- 
ing (shaking), the optical lattice potential. This pertur- 
bation couples vibrational levels in the lattice and creates 
an effective band structure through interference of tun- 
neling amplitudes for atoms in superposed vibrational 
states. The onset of instability is detected by directly 
observing the sudden growth of narrow momentum com- 
ponents of the condensate wavefunction at half the lattice 
recoil momentum, corresponding to perturbation modes 
which modulate the condensate wavefunction with a pe- 
riod of two lattice spacings. The critical onset threshold 
with shaking amplitude and frequency is measured and 
compared to the behavior expected from calculated band 
structure and numeric evolution of the GPE. 

To observe these effects, samples of 5 x 10^ bose- 
condensed atoms were prepared in the \F = 

2,rriF = 2) hyperfine state with no observable thermal 
coimonent by evaporative cooling in a triaxial TOP trap 
1^ [ 3 . After evaporation, the magnetic trap frequen- 
cies were relaxed to (ojqp,iiJnqp,oJz) = 27t( 55, 32, 28)Hz, 
and the condensate was adiabatically loaded into a set 
of 1,2 or 3 orthogonal standing waves derived from a 



low, 1.064/rm Nd:YAG laser injected by a stable seed 
laser (Lightwave NPRO model 125). Two standing waves 
made a 45° angle with the vertical, while the third 
was horizontal; each standing wave was formed by one 
retroreflected beam focused to 100/xm at the location 
of the condensate, and separated in frequency from the 
other standing waves by 80MHz. All beam intensities 
were actively stabilized to < 1% RMS over 2s. Adiabatic 
loading was achieved by linearly ramping the standing 
wave intensities over 400ms. To enable dynamic transla- 
tion of the lattice potential, an electrooptic phase modu- 
lator was placed in the retroreflection path of each stand- 
ing wave. In this work we used sinusoidal phase modu- 
lation = 4>ocos{uj£)t), which translated the lattice 
potential in the lab frame according to Viatt{x — xo(t)), 
where xo(t) = A(fo/47rcos(ojDt) is the modulation am- 
plitude and ljd its frequency. We calibrated the lat- 
tice by locating the heating resonance corresponding to 
phase modulation at the lattice vibrational frequency 
cuvib, where IkOyn, is the energy difference between the 
lowest two band centers. When multiple standing waves 
were used, their intensities were adjusted to equalize the 
vibration frequencies along all directions. 

After loading and interrogating condensates in the op- 
tical lattice, the magnetic trap and lattice were switched 
off in 200/is, and the atoms were absorbtively imaged @ 
after a 30ms time-of- flight. The condensate momentum 
distribution was inferred from the matter wave interfer- 
ence patterns of the released atoms. Without modula- 
tion, we observed interference peaks at momenta of nhk, 
where n is an integer, k is the inverse lattice vector, indi- 
cating phase-coherent localization in the lattice potential 
wells. Here, | k |= 27r/a = 2kught, where a is the lattice 
spacing and kught is the wavevector of the light forming 
the optical lattice. The relative peak heights were de- 
termined by the on-site wavefunction of the condensate, 
whose center of mass oscillated(Eig^) when the lattice 
was shaken near resonance. 

Modulating at a frequency below resonance in a one- 
dimensional lattice resulted in the growth of narrow 
peaks at momenta of ±fi.A:/2(Fig^), which also oscil- 
lated in phase with the drive (Eig^). The presence 
of peaks spaced by hk/2 indicates periodicity of the con- 
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FIG. 1: Absorption images of the atomic interference fol- 

lowing lattice modulation and time-of flight expansion, (a) 
Atoms are released at the turning-point of the lattice trans- 
lation, after a variable number of modulation cycles, Nc.yc- 
For short modulation times, peaks occur only at p = 0 and 
p = ±hk. For longer times, peaks develop at p = ±hk/2. (b) 
Release occurs within the 76th cycle, spaced by 7t/4 in phase; 
the center of mass oscillates nearly in phase with the applied 
drive, (c) Simultaneous doubling in a 2D lattice achieved by 
applying the same modulation as in (a,b) to two orthogonal 
standing waves, (d) The effect of modulation above reso- 
nance in a ID lattice; for small amplitudes, the p=0 peak 
develops a shoulder. At larger amplitudes, (e) broad peaks 
are apparent between p = 0 and p = ±hk. For visibility, (c) 
shows atomic absorbtion; (a,b,d,e) show atomic column densi- 
ties. Modulation and lattice parameters were {a,h,c)u>vib = Stt 
10.4kHz, xo = .04a, ujd = 2tv 7.6kHz, (d) uJvih = 2tv ll.OkHz, 
xq = .035o, lod = Stt 12.5kHz. (e) uJvib ~ Stt ll.OkHz, 

Xo = .037a, cjd = Stt 12.5kHz. 



densate wavefunction on a lengthscale of twice the lattice 
spacing. These peaks became visible only beyond a crit- 
ical frequency-dependent drive amplitude and following 

a sufficiently long modulation time. For a typical lattice 

1 2 

depth of 9Erec, where Erec = onset occurred after 
Sms of modulation with amplitude greater than 0.01a at 
a detuning of S = lod — tUmb = —2 tt 2 kHz. The ampli- 
tude for critical onset of this momentum class decreased 
with decreased detuning from resonance (Figl^J. Follow- 
ing onset, peaks were discernable under drive for ~ 800 
cycles (100ms), after which they were obscured by a loss 
of contrast in the interference pattern. We have found 
similar behavior at negative detunings for lattice depths 
from 3.9 to 19.5Erec- 

We have investigated the dependence of the critical 
drive parameters on the external trapping geometry by 
varying the envelope trapping frequencies by 30%. This 
caused no shift in the onset amplitude at a fixed detun- 
ing to 10%. We have also observed period-doubling in 
two-dimensional lattices of the same depth with one or 
both (Fig^) axes shaken. Though a decreased contrast 
in the interference patterns was observed, there was no 
measurable shift in the period-doubling onset amplitude 
at a fixed detuning. In three-dimensional lattices, loss of 
contrast precluded observation of coherent populations 
at IpI = hk/2. 

This sudden transition to non-zero momentum states 
is made possible by the effect of coupling the ground 
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FIG. 2: Critical drive strengths and detnnings. For 5 < 0, the 
period-doubling onset amplitnde decreases near resonance. 
For 5 > 0, the p=0 peak develops a shoulder, and the 
p = ±hk/2 classes are populated diffusely; these onsets are 
marked by open triangles. The inset shows relative peak 
atomic density at [p | = 0,hk/2 extracted from absorbtion 
images as a function of drive amplitude. Modulation pulse 
lengths for this data were 10ms. 



and first excited bands via near-resonant translation of 
the lattice potential. However, this perturbation acts 
equally on each lattice site and cannot itself induce mod- 
ulations of the wavefunction at twice the lattice spacing. 
Interparticle interaction, however, does allow for sponta- 
neous breaking of this discrete translational symmetry, 
and permits coupling between zero and non-zero quasi- 
momentum states. Accordingly, we consider as candidate 
mechanisms for this effect those involving the combined 
influence of interaction and lattice modulation. 

In the absence of modulation, Bloch states close to 
the first band center are stable against decay into other 
modes [^. Since the second band displays a compara- 
tively large, inverted tunneling rate, small modulation- 
induced admixtures of this state will dramatically alter 
the form of the single-particle dispersion, and the stabil- 
ity of Bloch waves must be reconsidered. For the case 
of far off-resonant modulation, band-mixin g ef fects have 
previously been described^ and observed|inj for ther- 
mal atoms. This effect becomes more apparent if we 
consider the two lowest bands in a ’dressed’ basis set 
similar to that in [^. Here, each basis state is a time- 
dependent superposition of delocalized ground and ex- 
cited states with the same quasimomentum, chosen to 
make the single-particle Hamiltonian time-independent 
(in a rotating wave approximation) and diagonal. 

The resulting band structure (Fig. |3 can be used to 
examine stability of stationary states in the presence of 
interaction. As modulation amplitude is increased, the 
lower dressed band is flattened by its growing admix- 
ture of the bare excited band, and eventually inverts its 
curvature at the band center. In analogy with results 
for lowest band Bloch wave stability in the tight-binding 
limit iiHH, one might postulate the emergence of 
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a dynamic instability of the q=0 state near the drive 
strength sufficient to invert the curvature of the lower 
dressed band. This criterion is equivalent to developing 
negative effective mass, and may be considered analogous 
to the case of attractive interaction with positive effec- 
tive mass. This prediction shows a qualitative agreement 
with the experimental data (Fig. for ^ < 0. 

To test this hypothesis, we performed linear stabil- 
ity analysis of the lattice-periodicity Bloch states for the 
lower dressed band in the tight-binding limit of the GPE, 
following methods in ii Adding a small perturbation 
of well-defined momentum to Bloch waves in this band, 
we looked for a parameter region where the spectrum 
of excitations includes imaginary values, signaling expo- 
nential growth or decay of perturbations to the wave- 
function. For negative detunings, one finds that zero 
momentum Bloch states first become dynamically unsta- 
ble when the modulation amplitude becomes sufficient 
to make the energy at the band edge equal to that at 
the band center Neglecting higher-order tunneling 
contributions, this condition is equivalent to demanding 
that the probability for an atom to hop to nearest neigh- 
bor sites is zero, while the hopping rate to next-nearest 
neighbor sites may remain nonzero. Close to onset, ex- 
ponential growth occurs only for perturbations with mo- 
mentum close to ±hk/2, consistent with experimental 
observations of sharp peaks at that momentum. The 
predicted threshold amplitudes for instability are smaller 
than those necessary to invert the effective mass, and 
show a better quantitative agreement with experimental 
data (Fig. at large negative detunings. This inter- 
pretation fails for small detunings, where the interaction 
energy exceeds the dressed band gap and introduces a 
non-negligible level mixing. 

We note that the curves shown in FigEI were cal- 
culated with knowledge only of the measured lattice 
vibration frequency, from which we calculate the un- 
perturbed band structure, tunneling rates for each 
band tr, and Wannier wavefunctions 4>a-{x). Coupling 
between bands is modeled by converting to a frame 
moving with the lattice, resulting in a single particle 
hamiltonian H = ^ -|- Viatt{x) + mxxo{t) |^. Here 
the effect of phase modulation is replaced by an iner- 
tial force, which gives an inter-band coupling strength 
hQ = muj'^xo f drd)ga;<i)i. The on-site interaction energy 
is estimated by approximating the on-site wavefunction 
transverse to the lattice ip(r±) by a Thomas-Fermi pro- 
file corresponding to measured transverse trapping fre- 
quencies and and an estimated average occupation num- 
ber. The interaction energy for atoms in bands a and 
a' is then given by Co-o-' = 

, where $(r(r) = tp{v i_)4>a{x), as is the s-wave scatter- 
ing length, and m the atomic mass. For simplicity, we 
approximate the on-site interaction energy Uo-a-> as inde- 
pendent of site and occupation number; more accurate 
treatments may be found in min 

To model the dynamics of the experiment, we have 
numerically evolved the CPE in the tight-binding limit 




FIG. 3: Renormalized band structure for negative (a,b) and 
positive (c) drive detuning <5, as calculated by diagonalizing 
the single particle hamiltonian in the dressed basis, includ- 
ing only the two lowest bands. Solid (dashed) lines indicate 
the energy of the renormalized ground (excited) band. These 
were calculated for a lattice depth of lO.SArec, and inter-band 
coupling strengths h^l of i)0, ii)0.3 and iii)0.9 Emc- Detuning 
is labeled in units of Erec/k. 

for two coupled bands. Using the approach from [ihlll^ . 
we describe the condensate wave function with a complex 
site-dependent amplitude Cj^a- for each (bare) vibrational 
level a, (r|fl>) = J2j ffi>d>o-(r — ja). We substitute (r|fl>) 
into the CPE and numerically evolve the amplitudes Cj^a 
according to 

Cj — Uj Cycr H ^ 

-f ^ ^ -f hEl COS {uj ^ 

( 1 ) 

where 

U/ = - xo{t) f + hwyibSi^^ + mxo{t)ja (2) 

represents the external trapping potential, band energy 
and the intra-band effect of the inertial force. We initial- 
ize the calculation with a lower band wavefunction solved 
to be stationary in absence of modulation. The chemical 
potential is chosen to reproduce the number of atoms N 
in the experiment. The calculation was performed with 
an external potential comparable to that used in the ex- 
periment, and number of sites chosen to safely exceed the 
number occupied in the experiment. 

The parameters for onset of the ±hk/2 peaks at neg- 
ative drive detunings show good agreement with those 
observed in the experiment (Fig. 2J, and agree well with 
the prediction for instability at the center of the low- 
est dressed band described above for large detuning. We 
have performed the same numerical evolution for an ex- 
cited state with a positive tunneling amplitude; in this 
case, no period-doubling was observable. 

To investigate the importance of the external trapping 
potential, we compare period-doubling onsets with and 
without LOext- To reduce numeric artifacts stemming from 
the large inertial force and boundary effects, this compar- 
ison was done in the stationary frame, to first order in 
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FIG. 4: Comparison of experimental and numeric data for 
drive parameters at onset. Boxes show experimental data 
points for onset under same conditions as FigH Solid line 
shows analytic result for equal single-particle energy at band 
center and edge, which occurs for zero drive amplitude when 
bands cross. Dashed (dotted) line shows line of infinite effec- 
tive mass at band center (edge). Between dashed and dot- 
ted lines effective mass is positive both at band center and 
edge. Circles show onset in numeric evolution of the CPE for 
Jo = 0.024, Ji = -0.28, Uoo = Uw = (4/3)t/n = 2 x lO"'^ in 
units of Erec\ u!ext = 2n X 45Hz and N = 5 x 10^. Note that 
the numeric data is a zero-parameter comparison to experi- 
mental data generated from the measured uj^ib- 

the lattice displacement amplitude xq. Without the ex- 
ternal potential, periodic boundary conditions were used 
with a site number comparable to the experiment (40). 
In this case, we find a small initial seed population at 
p = ±hk/2 grows rapidly once a sufficient modulation 
strength is applied. Without a seed population, observ- 
able peaks at ±hk/2 do not occur for experimentally rel- 
evant timescales. With an external trapping potential 
this seed population is guaranteed by the finite extent 
of the condensate wavefunction, and numerically results 
in an observable population at p = ±hk/2 in several ms, 
consistent with experimental observations. This suggests 
that, while the mechanism for period-doubling exists on 
a lattice with discrete translational symmetry, the on- 



set of the period-doubling population is hastened by an 
external potential through its effect on the initial state. 

The analysis above applies to modulation performed 
below resonance, where the influence of the avoided cross- 
ing in Figiniis felt most strongly at the band edge. Exper- 
imentally, modulation performed at frequencies higher 
than the vibrational resonance first results in the ap- 
pearance of a shoulder around the p=0 peak, followed 
by growth of broad peaks around p = ±Hkj2 (Fig. Eije). 
The critical drive amplitude necessary for this behavior 
is higher than for the period-doubling onset at an equal 
detuning of opposite sign (Fig. [21 . While such behavior 
might be expected from the band structure in curve ii of 
FigEt, the details of its onset and dynamics are not yet 
fully understood. 

In conclusion, we have experimentally observed a 
period-doubling instability of a condensate held in 
a shaken optical lattice potential, and determined the 
threshold modulation conditions necessary for the on- 
set of this behavior. A simple physical interpretation 
provides qualitative agreement with the onset of the ob- 
served dynamic instability, while a numerical simulation 
of the effect based on the GPE is in quantitative agree- 
ment with this picture. 

A new feature in this observation of dynamic insta- 
bility is the controlled introduction of instability at zero 
quasimomentum and the ability to directly observe the 
growth of a well-defined unstable mode. Due to the depth 
and dimensionality of the lattices used, it is interesting 
to consider the effect of fluctuations of the long range co- 
herence of the condensate near the onset of this instabil- 
ity. In the dressed-state picture, the modulation-induced 
flattening of the lower band at q=0 further increases the 
already large ratio of local interaction energy to tunnel- 
ing amplitude, pointing to an increasing importance of 
fluctuations and limited validity of a mean field descrip- 
tion. In this light, a description beyond that of the GPE 
may warrant future study. 
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